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Abstract. We discuss the derivation of the analytic properties of the cross- 
power spectrum estimator from multi-detector CMB anisotropy maps. The 
method is computationally convenient and it provides unbiased estimates under 
very broad assumptions. We also propose a new procedure for testing for 
the presence of residual bias due to inappropriate noise subtraction in pseudo- 
C'i estimates. We derive the analytic behavior of this procedure under the 
null hypothesis, and use Monte Carlo simulations to investigate its efficiency 
properties, which appear very promising. For instance, for full sky maps with 
isotropic white noise, the test is able to identify an error of 1% on the noise 
amplitude estimate. 



PACS numbers: 95.75Pq, 98.80Es, 02.50Ng, 02.50Tt 

1. Introduction 

The cosmic microwave background (CMB) provides one of the most powerful ways 
of investigating the physics of the early Universe. The main CMB observable is the 
angular power spectrum of temperature anisotropy, which encodes a large amount of 
cosmological information. In the last decade, important advances in the measurement 
of the CMB angular power spectrum took place; this resulted in relevant progress 
in our understanding of physical cosmology. CMB temperature anisotropics were 
first detected by the COBE satellite in 1992 [J. This discovery fuelled a period 
of intensive experimental activity, focused on measuring the CMB power spectrum 
on a large range of angular scales. A major breakthrough was made in the past 
few years, when the MAXIMA and BOOMERanG balloon-borne experiments 
independently produced the first high-resolution maps of the CMB, allowing a clear 
measurement of a peak in the power spectrum, as expected from theoretical models 
and previously detected by the ground based experiment TOCO 0]. Since then, 
many other experiments have confirmed and improved on these results: DASI0, 
BOOMERanG-B98 01313, BOOMERanG-B03 |ni[II3IIll|T2[l3|, VSA [H], Archeops 
[El, CBI [ISj, ACBAR[n], BEAST J^. Most notably, the NASA satellite mission 
WMAP, whose first year data were released in February 2003 (JHI and references 
therein), provided the first high-resolution, full sky, multi-frequency CMB maps, and 
a determination of the angular power spectrum with unprecedented accuracy on a 
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large range of angular scales. Much larger and more accurate data sets are expected 
in the years to come from ESA's Planck satellite. 

In this paper, we shall concentrate on extracting the CMB power spectrum from 
full sky maps with foregrounds removed. We shall focus mainly on techniques for 
dealing with noise subtraction. In principle, and for Gaussian maps, noise subtraction 
can be performed by implementing maximum likelihood estimates. It is well known, 
through 20, HI, that maximum likelihood estimates require for their implementations 
a number of operations that scales as Np^^, Npix denoting the number of pixels in the 
map. For current experiments, Npix ranges from several hundred thousands to a 
few millions, and thus the implementation of these procedures is beyond computer 
power for the near future. Many different methods have been proposed for producing 
computationally feasible estimates; here we just mention a few of them, and we refer 
the reader to for a more complete discussion on their merits. Some authors 
have introduced special assumptions on the noise properties and symmetry of the 
sky coverage, to make likelihood estimates feasible; see, for instance, [231 1^ ESI ESI ■ 
Reference |2Z| adopted an entirely different strategy, extracting the power spectrum 
from the 2-point correlation function of the map. Others have used estimators based on 
pseudo-C; statistics and Monte Carlo techniques [2111221 j or based on Gabor transforms 
|30| . For multi-detector experiments, an elegant method, based on spectral matching 
to estimate jointly the angular power spectrum of the signal and of the noise, was 
proposed in [31] . Pseudo-Cf estimators were adopted by the WMAP team [32| , which 
used the cross-power spectrum estimator and discussed the best combination of the 
cross-power spectrum obtained from single couples of receivers. 

Our purpose in this paper is to derive some analytic results on the cross-power 
spectrum estimator, to perform a comparison with standard pseudo-C^ estimators, 
and to propose some testing procedures on the assumption that any noise bias has 
been properly removed, which is clearly a crucial step in any estimation approach. 
We shall also present some Monte Carlo evidence on the performance of the methods 
that we advocate. The plan of this paper is as follows. In Section [21 we derive the 
analytic properties for the cross-power spectrum estimator and we compare them 
with equivalent results on standard pseudo-C^ estimators. In Section [31 we propose a 
procedure (the Hausman test) for verifying appropriate noise subtraction in pseudo- 
Cf, estimators, and we derive its analytic properties. In Section [3 we validate our 
results by using Monte Carlo simulations, which are also used to test the power of our 
procedure in the presence of noise which has not been completely removed. In Section 
ISlwe review our results and discuss directions for future research. 

2. Power spectrum estimators 

The CMB temperature fluctuations ■^(6', (/>) can be decomposed into spherical 
harmonic coefficients 



Jn J 

If the CMB fluctuations are Gaussian distributed and statistically homogeneug, 
as suggested by the latest experimental results (see for instance [331 [^[3S]), then 
each atra is an independent Gaussian complex variable with 




(1) 



{aim) = 



(2) 
(3) 
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and all the statistical information is contained in the power spectrum Ci. 

In the following we describe two procedures for estimating the CMB angular 
power spectrum: the standard pseudo-C^ estimator, sometimes labelled the auto- 
power spectrum j32l . and the cross-power spectrum. As a first step, we shall assume 
handling of full sky maps with isotropic, not necessarily white noise. 



2.1. The standard pseudo-Ce estimator 

Pseudo-Cf estimators are very useful in computing the power spectrum because they 
are fast enough to be used on large data sets such as WMAP and Planck. The standard 
pseudo-Cf estimator has been thoroughly investigated in the literature, taking also into 
account some important features of realistic experiments such as partial sky coverage 
and systematic effects 2 8^ . The starting point is the raw pseudo-power spectrum 
defined as 

1 ^ 

= 2e+l ^ di„^dem (4) 

m——£ 

where dim are the spherical harmonic coefficients of the map. 

In the absence of noise and for a full sky CMB map, dgm — and is an 
unbiased estimator of Ci (the angular power spectrum of the signal) with mean equal 
to Cg and variance equal to 2Ce/2£ + 1; also, {2£ + \)Cf' /C^ is a x^— distributed 
variable with v = 2£ -\- 1 degrees of freedom. 

In the presence of noise, it is not difficult to see that this estimator is biased. If 
we assume, as usual, that noise is independent from the signal, we have 

and 

(df ) = Q + . (6) 

Now the common assumption is to take as determined a priori, for instance by 
Monte Carlo simulations and measurements of the properties of the detectors; we shall 
discuss later how to test the validity of this assumption and/or make it weaker. Under 
these circumstances, the power spectrum estimator is naturally defined as 

Cl = CP - = 21+1 ^ dirndlm - . (7) 

m——£ 

Of course, if the estimate of the noise power spectrum is not correct, the estimator 
will be biased. For a multi-channel experiment, we generalize equationElby averaging 
the maps from each detector and then computing the power spectrum of the resulting 
map. A more sophisticated approach would be to use weighted averages, with weights 
inversely proportional to the variance of each detector, but we shall not pursue this 
idea for the sake of brevity. In view of equations and in the presence of k channels 
with uncorrelated noises we can write 



Unbiased estimation of an angular power spectrum 



4 



where i is the detector index and afj^ are the noise spherical harmonics coefficients. 
Assuming that our noise estimation is correct, we obtain for the expected vahie and 
the variance 



< Ce >= Co 



and 



21- 



i=i i=i 



(9) 



(10) 



It should be noted that in equation 1101 the value of Cf is taken as fixed, and in this 
sense we are underestimating the variance by neglecting the additional uncertainty 
due to the estimation of the noise properties. 



2.2. The cross-power spectrum 

The pseudo-Cf estimator presented in the previous subsection is computationally very 
fast and simple to use, but it is prone to bias if noise has not been appropriately 
removed. It is thus natural to look for more robust alternatives, yielding unbiased 
estimates even in the presence of noise with an unknown angular power spectrum. 
For this purpose, we now focus on the cross-power spectrum, which is defined, for any 
given couple of channels (i, j), as 



CI- 



21 



-J] 



It iss easy to show that 



and 



Var 



{cf} 



21+1 



Ci 



(11) 



(12) 



(13) 



For the details of the calculations see the appendix. Let us now consider the most 
general case with k detectors; this means that we can construct k{k — l)/2 different 
couples of channels. For each of them we can calculate the cross-power spectrum and 
then take the average; thus the cross power spectrum becomes 

fc-l k 



Cf, = 



fc(fc- 1) 



i=l j=i+l 



Again, the resulting estimator is clearly unbiased, {Ct 

by 



(14) 

Cf. Its covariance is given 



Var 



fe-1 



Var 



k{k 
4 



T)E E ^- 

' i=l j=i+l 



k-1 



fc2(fc- 1)2 
4 



E E ^-{^;^} 



(15) 



i—1 j—i^l 



k^{k-i) 



-2 [2Cov {cf , } + 2Cov {cf , Cf } 
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In order to evaluate this quantity, the first step is to consider the covariances among 
different pairs («, j), {h, k). For k channels we can construct c = k{k — l)/2 different 
couples and t — c(c — l)/2 covariance terms, which are 

"^—Cf for h ^ i,j and k ^ i,j 



Cov 



{•-t^^i I - I ^ |c| + iQCf } for = z or j and k^i.j . 



The next step is to consider how many times we have the term, for 

each i = 1, ...,/c. This term appears when one of the two index of a couple is equal 
to one of the two index of another couple. This leaves (fc — 1) possible values for the 
second index in the first couple, and (k — 2) possible values for the second index in 
the second couple; finally we have a factor 1/2 to take into account symmetries, that 

is, the fact that Cov^C\\C'l^^ = Cov ^C^^ ,0^'^ (equivalently, we could drop the 

factor 2 which multiplies the covariance terms in eauation llBf) . The result is that the 
single term CfC^'/2 appears (fc — l)(fc — 2)/2 times. 

By plugging in equation 1161 into equation 1161 we obtain 



Var 



{^'] = ^? - > t + E t CtcP \ . (17) 



^ ^ i=l j=i+l 



It can be verified that for fc = 2, equation ^| reduces to equation E| It is 
interesting to compare this result with the variance of the classic pseudo-C^ estimator. 
We can write immediately 



Vc 



i—1 i—1 J— 2+1 



Considering the case where ' — for all the channels, we obtain 

y^r{a}^Var{c,}^^^[^^^iCn^] . (19) 

Hence, if noise has the same power spectrum over all channels, then the standard 
estimator is always more efficient, although clearly the difference between the two 
estimators becomes asymptotically negligible as the number of detectors grows (it 
scales as fc~^). 

We have thus shown that the cross-power spectrum estimator provides a robust 
alternative to the classical pseudo-C^ procedure, in that it does not require any a 
priori knowledge of the noise power spectrum. We shall argue that cross-power 
spectrum estimates can be extremely useful even if different procedures are undertaken 
to estimate the angular power spectrum; indeed, in the next section we discuss how 
to test the assumption that noise has been appropriately removed from the data from 
a multi-channel experiment. 

3. The Hausman test 

In the previous section, we compared the relative efficiency of the two estimators 
C£,Ci, in the case where the bias term in had been effectively removed. In this 
section we propose a testing procedure to verify the latter assumption. Consider the 
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random variable = — ; if is unbiased, then it is immediate that Ge has 
mean zero, with variance 

Var ^Ci - q| = Var i^Ci^ + Var |c^| - 2Cov jc^, C^j , (20) 

where 

k — l k 

Cov{Ce,C,}^—^Yl E Cov{c,,Cl^} . (21) 

^ ' i=l j=i+l 

In the appendix we show that, for a single couple we have 

= + } . (22) 

Now we use equation 1221 in eg uation 12 II and we obtain 

k ^ k—\k 



^ ^ { C| + ^ E + E E_ ) . (23) 

Therefore 



fc3(fc- 1) 
i=l ^ ^ 1=1 j=i+l 



k^(k ^, 
i—1 2=1 



The special case = ... = Cf'' gives 

2 1 
2^+1 F(fc- 1) 
Thus, for a fixed £ we can suggest the statistic 

^1/2 



1/2 

2, 



(k ^, 

^ ' i=l j=i+l 



(26) 



as a feasible test for the presence of bias in Ce- By a standard central limit theorem, 
we obtain that 

Hi -^'^ iV(0, 1) as £ ^ cx) (27) 

where -^'^ denotes convergence in distribution and A^(0, 1) represents a standard 
Gaussian random variable. In words, for reasonably large £ the distribution of Hi 
is very well approximated by a Gaussian, provided that Cg is actually unbiased; on 
the other hand, if this is not the case the expected value of Hi will be non-zero. This 
observation suggests many possible tests for bias, using for instance the chi-square 
statistic Hi (a value of Hf larger than 3.84, the chi-square quantile at 95%, would 
suggest that bias has not been removed at that confidence level). In practice, however, 
we have to focus on many different multipoles, 11=1, ■■■,L, where L depends on the 
resolution of the experiment and its signal to noise properties. It is clearly not enough 
to consider the whole sequence {Hi}^^^ ^ l ' check for the values above the 
threshold, as this does no longer correspond to the 95% confidence level (it is obvious 
that, if P(iJ| > 3.84) ~ 5%, then P{m&M=i,...,L Hj > 3.84) >> 5%, the exact value 
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being difficult to determine). To combine the information over different multipoles 
into a single statistic in a rigorous manner, we suggest the process 

-, [Lr] 

where [.] denotes integer part. Of course, other related proposals could be considered; 
for instance we might focus on weighted versions of i?L(r), to highlight the contribution 
from low multipoles, where it is well known that there are problems with non-maximum 
likelihood estimators. This modification, however, would not alter the substance of 
the discussion that follows. 

We note first that Bi^ir) has mean zero; indeed, 

^ [Lr] 

<Si(r)>=^^<i/, >=0. (29) 

Also, for any ri < r2, as L — s- cx), 

[in] [Lt2] [in] 
Cov{BL{ri),BL{r2)} ^ J J2 J2 < ^^^^ >= jY. < >^ ""^ ' (3^) 

As r varies in [0, 1], B^ir) can be viewed as a random function, for which a functional 
central limit theorem holds; in fact, because BL{r) has independent increments and 
finite moments of all order, it is not difficult to show that, as L ^ oo, 

BLir)^Bir) (31) 

where =4» denotes convergence in distribution in a functional sense (see for instance 
|35)1: this ensures, for instance, that the distribution of functionals of B]^{r) will 
converge to the distribution of the same functional, evaluated on B{r). Also, B{r) 
denotes the well known standard Brownian motion process, whose properties are 
widely studied and well known: it is a Gaussian, zero-mean continuous process, with 
independent increments such that 

B{r2)-B{ri)^N{0,r2-ri) . (32) 

In view of eauation l31l and standard properties of Brownian motion, we are for instance 
able to conclude that 

lim P I sup BL{r) > x \ = P I sup B{r) > x \ ^ 2P{Z > x), (33) 

L^oo [re[0,l] J [re[0.1] J 

Z denoting a standard (zero mean, unit-variance) Gaussian variable (see for instance 
|37|V This means that to determine approximate threshold values for the maximum 
value of the sum X^fcTi ^ varies between zero and one, the tables of a 

standard Gaussian variate are sufficient. Likewise, the asymptotic distribution of 
P{sup^ \BLir)\ > x} is given by 

lim P \swp\BLir)\ > x\ = (34) 



E 



k— — oo 



(z + Akx)"^ \ f (z + 2x + Akxf 

exp I — — - cxp ' 



dz 



Monte Carlo simulations have confirmed that equation 1331 and equation 1341 provide 
accurate approximations of the finite sample distributions, for L in the order of 10'^. 
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4. Effect of noise correlation 



In order to consider the effect of correlated noise we start discussing tlie simplest case 
with two detectors. The presence of correlated noise can be inserted by rewriting 
equation |S1 as 

4m = + ftfo,! + C^lrn (35) 

where c'g^ is independent from aim-, and a^. Under these circumstances, it is 
clear that both Cg and Gg will be biased; however, their difference G^, used in the 

(12) 

Hausman test, is not affected at all due to cancellations of all the terms involving c^^ : 



m— — ^ 

i('lo^i|2 + |„^2|2 , r, /Vi-Af2N , / Wl JV2N-(12) |J12)|2 

(\ |2 , /-iVi , -JVa , o-(12)n , JVi-TVa , / iVi , N^\-{\1) , i (12) ,2^ 

(^la^ml + a£m(%^ + + ^c\J) + a^„>^^ + + ai,^l)c)J + |c)„/ 1 J j 

In the more general case with k detectors with noise correlations which varies 
from pair to pair, this is no longer true. In fact, by completely analogous arguments, 
it can be shown that some extra terms involving cross-products of the form a^c^^ 
will remain in Gi. These terms, however, have zero expected value, and thus will 
affect only the variance of H^. In other words, the previous approach can go through 
unaltered, provided we have available a reliable estimate of the variance of Gi. These 
issues are investigated by means of Monte Carlo simulations in the next section. 



5. Monte Carlo simulations 



To verify the validity of the previous analytic arguments, we present in this section 
some Monte Carlo simulations. As a first step, we generate some Gaussian, full sky 
CMB maps from a parent distribution with a given power spectrum, which corresponds 
to a standard KG DM model with running spectral index; the values of the parameters 
are provided by the WMAP best fit, that is f^b/i^ = 0.02262, Q.cDMh? = 0.10861, 
n{k = 0.05Afpc-i) = 1.04173, exp{-2T) ^ 0.69879, dn/dlnk = -0.01618, ampik = 
O.OSMpc"^) = 0.86746, h = 0.73070. In order to include the effect of a finite resolution 
of the detectors, we simulate the maps using a beam of 12' FWHM. Then we considered 
two channels and added random Gaussian noises realizations to each of them; noise 
is assumed to be white and isotropic with RMS amplitude per 7' pixels of 55^K and 
GbfiK respectively for the two channels. The input power spectra used are shown if 
figure n We start considering full sky maps. From each CMB realization we compute 
both the cross-power spectrum and the auto-power spectrum, for I = 2, ...,L = 1300. 
We generated 1000 maps, and we start by presenting the Monte Carlo values for the 
variances of the cross-spectrum and auto-power spectrum estimators, together with 
the variance of their differences. Results are shown in figures EBI they are clearly in 
extremely good agreement with the values that were obtained analytically. 
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Figure 1. Input power spectra used in the simulation. The RMS of the noise 
per 7' pixel is 55 and Qb/iK for channel 1 and channel 2 respectively. 



auto — power spectrum standard deviation 
3000 1 I 
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Figure 2. Standard deviation of the auto-power spectrum estimator. In grey 
we show the results obtained from 1000 Monte Carlo simulations, while the black 
line is obtained from equation 1171 
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Figure 3. Standard deviation of the cross-power spectrum estimator. In grey 
we sliow the results obtained from 1000 Monte Carlo simulations, while the black 
line is obtained from equation 1101 




Figure 4. Standard deviation of the difference between the auto- and the cross- 
power spectrum estimator. In grey we show the results obtained from 1000 Monte 
Carlo simulations, while the black line is obtained from equation 1241 
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Table 1. Threshold values under the null hypothesis 







S2 


S3 


68% 


0.919 


1.533 


0.596 


95% 


1.818 


2.550 


1.926 


99% 


2.401 


3.308 


3.842 



We now focus more directly on the efficiency of the Hausman test in identifying 
a residual bias in the auto-power spectrum. In order to achieve this goal, we simulate 
300 further maps with a noise power spectrum , and we compute the auto-power 
spectrum using a modified version of equation [7| 

^ 1 ^ 

Ce = 2g ^ I '^irndlm — fnCf ■ (37) 

rn——i 

In this way we simulate a wrong estimation of the noise power spectrum. 

Then, for a fixed /„, we compute and B^^r) for each simulation. We consider 
the three test statistics si — sup^i3j;,(r), S2 — sup^\Bi^{r)\ and S3 = Bj^{r)dr, 
and the threshold values for the 68%, 95% and 99% probability. We used a thousand 
independent simulations with the value /„ — 1,, corresponding to the case where our 
a priori knowledge of noise is correct, to tabulate the empirical distributions under 
this null hypothesis; results are reported in tabled 

We then go on to compute si,S2 and S3 under the alternatives /„ =/= 1; the 
percentages of rejections provide an estimate of the power of these procedures in 
detecting a bias. Results are reported in tables I^BI and are clearly very encouraging: 
the S2 and S3 test statistics enjoy 100% power even in the presence of a mere 1% 
misspecification of the noise angular power spectrum. Note that, as expected, si is a 
unidirectional test, that is, it has no power in the case where noise is overestimated 
(/n > 1); however, for such circumstances it would suffice to consider s[ = inf^ BL{r) 
to obtain satisfactory power properties. In general, S2 and S3 should clearly be 
preferred for their robustness against a wider class of departures from the null. 



Table 2. The power of the test si 



u 


0.99 


0.995 


0.998 


0.999 


1.001 


1.002 


1.005 


1.01 


68% 


1.00 


1.00 


0.79 


0.54 


0.15 


0.08 


0.01 


0.00 


95% 


1.00 


0.99 


0.39 


0.14 


0.01 


0.00 


0.00 


0.00 


99% 


1.00 


0.96 


0.16 


0.05 


0.00 


0.00 


0.00 


0.00 


Table 3. The power 


of the test S2 










/n 


0.99 


0.995 


0.998 


0.999 


1.001 


1.002 


1.005 


1.01 


68% 


1.00 


1.00 


0.54 


0.27 


0.53 


0.82 


1.00 


1.00 


95% 


1.00 


0.95 


0.12 


0.04 


0.11 


0.43 


0.98 


1.00 


99% 


1.00 


0.76 


0.04 


0.00 


0.03 


0.12 


0.95 


1.00 
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Table 4. The power of the test S3 



fn 


0.99 


0.995 


0.998 


0.999 


1.001 


1.002 


1.005 


1.01 


68% 


1.00 


0.99 


0.57 


0.33 


0.45 


0.74 


1.00 


1.00 


95% 


1.00 


0.86 


0.15 


0.06 


0.11 


0.30 


0.95 


1.00 


99% 


1.00 


0.58 


0.03 


0.01 


0.02 


0.08 


0.76 


1.00 



5.1. Effect of partial sky coverage 

In order to study the effect of partial sky coverage on the Hausman test, we repeated 
the Monte Carlo analysis considering the patch observed by BOOMERanG, covering 
^ 2% of the full sky j3Slin|- We expect this to be a good limiting case, where any 
failures of the test due to partial sky coverage should clearly show up. 

The main effect of partial sky coverage is, as well known, to produce correlations 
among spherical harmonic coefficients, that can be interpreted as a reduction of the 
effective number of degrees of freedom in the power spectrum !28J . 

Results are reported in tables I5I7I We stress that the power of the Hausman 
test, although reduced (as expected), is still very satisfactory. For instance, a 
misspecification of the noise level of the order of 5% is detected 100% of the times by 
S2 and 99% by S3. 



Table 5. The power of the test si in the presence of partial sky coverage. 



fn 0.95 


0.96 


0.97 


0.98 


0.99 


1.01 


1.02 


1.03 


1.04 


1.05 


68% 1.00 


1.00 


0.99 


0.94 


0.70 


0.09 


0.03 


0.01 


0.00 


0.00 


95% 1.00 


1.00 


0.97 


0.75 


0.35 


0.01 


0.00 


0.00 


0.00 


0.00 


99% 1.00 


0.99 


0.91 


0.63 


0.17 


0.01 


0.00 


0.00 


0.00 


0.00 



Table 6. The power of the test S2 in the presence of partial sky coverage. 



/„ 0.95 


0.96 


0.97 


0.98 


0.99 


1.01 


1.02 


1.03 


1.04 


1.05 


68% 1.00 


1.00 


0.99 


0.88 


0.58 


0.42 


0.82 


0.98 


1.00 


1.00 


95% 1.00 


1.00 


0.95 


0.67 


0.26 


0.18 


0.54 


0.91 


0.99 


1.00 


99% 1.00 


0.99 


0.85 


0.50 


0.11 


0.08 


0.34 


0.78 


0.97 


1.00 



Table 7. The power of the test S3 in the presence of partial sky coverage. 



fn 


0.95 


0.96 


0.97 


0.98 


0.99 


1.01 


1.02 


1.03 


1.04 


1.05 


68% 


1.00 


1.00 


0.96 


0.80 


0.51 


0.37 


0.70 


0.94 


0.99 


1.00 


95% 


1.00 


0.96 


0.82 


0.53 


0.16 


0.13 


0.35 


0.73 


0.95 


0.99 


99% 


0.99 


0.92 


0.70 


0.33 


0.09 


0.05 


0.23 


0.55 


0.88 


0.99 
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5.2. Polarization and 1 /J noise correlated among different detectors 

We move forward, analysing a more realistic case including polarization measurements 
in the presence of 1/f noise correlated among different detectors. This is achieved 
generating time ordered data with a scanning strategy and detectors noise properties 
similar to those of BOOMERanG-B03, where correlations of the order of 10% are 
present (see table 7 in 0), and the 1/f noise knee frequency is ~ Q.QTHz (see figure 
21 in [H]). The sky maps are then obtained using the ROMA IGLS polarization map- 
making code 

Polarization measurements provide six power spectra that can be used separately 
or combined to obtain a more efficient detection of the noise bias. The optimal 
combination of polarization power spectra is under investigation and will be addressed 
in a future paper. Here, in order to illustrate the method, we simply average the Bi^{r) 
obtained from each power spectrum. 

Results are reported in tables ISITUl Once more the power of the Hausman test is 
reduced with respect to the full sky uncorrelated noise case, but is still satisfactory. 
For instance, a misspecification of the noise level of the order of 15% is detected 
100% of the times with 95% significance. 



Table 8. The power of the test si for polarization measurements with 1/f noise 
correlated among different detectors. 



/n 


0.80 


0.85 


0.90 


0.95 


0.98 


1.02 


1.05 


1.10 


1.15 


1.20 


68% 


1.00 


1.00 


1.00 


0.93 


0.79 


0.48 


0.29 


0.12 


0.05 


0.00 


95% 


1.00 


1.00 


0.89 


0.47 


0.21 


0.06 


0.04 


0.02 


0.01 


0.00 


99% 


0.99 


0.88 


0.45 


0.11 


0.05 


0.02 


0.01 


0.00 


0.00 


0.00 



Table 9. The power of the test S2 for polarization measurements with 1/f noise 
correlated among different detectors. 



/n 


0.80 


0.85 


0.90 


0.95 


0.98 


1.02 


1.05 


1.10 


1.15 


1.20 


68% 


1.00 


1.00 


0.99 


0.83 


0.65 


0.76 


0.90 


1.00 


1.00 


1.00 


95% 


1.00 


0.99 


0.73 


0.29 


0.12 


0.17 


0.39 


0.86 


0.99 


1.00 


99% 


0.99 


0.83 


0.38 


0.07 


0.04 


0.04 


0.07 


0.51 


0.93 


1.00 



Table 10. The power of the test S3 for polarization measurements with 1/f noise 
correlated among different detectors. 



fn 


0.80 


0.85 


0.90 


0.95 


0.98 


1.02 


1.05 


1.10 


1.15 


1.20 


68% 


1.00 


1.00 


0.98 


0.84 


0.72 


0.75 


0.89 


1.00 


1.00 


1.00 


95% 


1.00 


0.95 


0.61 


0.21 


0.14 


0.17 


0.36 


0.76 


0.96 


1.00 


99% 


0.97 


0.70 


0.25 


0.09 


0.04 


0.05 


0.07 


0.36 


0.82 


0.99 



The Planck experiment will provide full sky polarization maps with pixel 
sensitivity similar to that of BOOMERanG-B03. Such a wide sky coverage will allow 
us to reach unprecedented accuracy in the estimated power spectra. The application 
of the Hausman test to Planck simulated maps will be the subject of a forthcoming 
paper. 
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6. Conclusions 



We have discussed the analytic properties of the cross-power spectrum as an 
estimator of the angular power spectrum of the CMB anisotropics. The method is 
computationally convenient for very large data sets as those provided by WMAP or 
Planck and it provides unbiased estimates under very broad assumptions (basically, 
that noise is uncorrelated along different channels). It thus provides a robust 
alternative, where noise estimation and subtraction are not required. We also propose 
a new procedure for testing for the presence of residual bias due to inappropriate 
noise subtraction in pseudo-Cf estimates (the Hausman test). The test compares the 
auto- and cross-power spectrum estimators under the null hypothesis. In the case of 
failure, the more robust cross-power spectrum should be preferred, while in the case 
of success both estimators could be used, and the choice should result from a trade-off 
between efficiency and robustness. We derive the analytic behaviour of this procedure 
under the null hypothesis, and use Monte Carlo simulations to investigate its power 
properties, which appear extremely promising. We leave for future research some 
further improvements of this approach, in particular, the use of bootstrap/resampling 
methods to make even the determination of confidence intervals independent from 
noise estimation. Finally, the optimal combination of polarization power spectra 
and the application of the Hausman test to the Planck satellite are currently under 
investigation. 
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Appendix 

We recall the definition of the cross spectrum estimator (equation ^J: 

m—~i m— — £ 

1 ^ 

\ " f/r,! |2 , Ni- I -Ni , -Nj , _ Nj , Ni-Nj , _Ar. Nj~.\ 

= 27TT ^ + a-i^aern + a^^^^aim + aima^^ + ae^a^^ + a^^a^,^ + a^^a^^) j 

m— 1 

It is easy to see that all summands in equation lA.ll are uncorrelated (albeit not 
independent), with variances given by 

Var {2\a,„,,\^} ^ ACi,m ^ I, ...,£, Var {\a,o\''} ^ 2Ci (A.2) 
^^a^^air,^^ =Far|a^a£™| ^ C^'Ci,m ^ 0, ...,£ (A. 3) 

V^r {af^af^} = Var {af^^af^} = ^ , m = 0, £ (A.4) 
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whence we obtain 

^ ' rn=l 

We show now that, for a single couple (i, j), we have 
Indeed, 

Cov\dt,Cl^^ = 

= Gov ^{\aimf + atmatm + o-imaim + \aim?] ,{\aim? + at-mafm + a^^aft" + 



^ Cot; 



(2^+1)2 



where 



Now, 
Gov 



1 

^ ' m=l 

k 
i=l 



and likewise 

Gov [2 {\aemf + aemaZn + o^emaem + > 



4 <! C| + ^(Cf • + Cf ) + Icf -C;^ ^ . (A.IO) 



Hence, 



m— 1 ^ 



'(2£+l)2^ 

^ [cl + ^(C^f + + ^,C-^Cf } (A.11) 



as claimed. 



Unbiased estimation of an angular power spectrum 



16 



References 

Smoot GF et al 1992 ^pJ396 LI 
Hanany S et al 2000 ApJ 545 L5 
de Bernardis Pcet al 2000 Nature 404 955 
Miller A D et al 1999 ApJ 524 LI 
Halverson K W et al 2002 ApJ 568 38 
Netterfield C B et al 2002 ApJ 571 604 
de Bernardis P et al 2002 ApJ 564 559 
Ruhl J et al 2003 ApJ 599 786 



Masi S et al 2005 Preprint astro-ph/0507509" 
Jones W C et al 2005 Preprint astro-ph/0507494] 
Piacentini F et al 2005 Preprint astro-ph/0507507' 
Montroy T E et al 2005 Preprint astro-ph/0507514 
MacTavish C 3 et al 2005 Preprint astro-ph/0507503] 
Grainge K et al 2003 MNRAS 341 L23 
Benoit A et al 2003 ASA 399 L19 
Pearson T J et ai 2003 /IpJ 591 556 
Kuo CL et al 2004 ApJ 600 32 
O'Dwyer IJ et a/ 2005 ApJS 158 93 
Bennett C h et al 2003 ApJS 148 1 

Bond J R , Jaffe A H and Knox L E 1998 Phys. Rev. D 57 2117 

Borrill J 1999 AlP Conf Proc 476 277B 

Efstathiou G 2003 Preprint astro-pli/03075T5l 

Oh S P , Spergel D and Hinshaw G 1999 ApJ 510 550 

Wandelt B, Hivon E and Gorski K M 2001 Phys. Rev. D 64h3003W 

Wandelt B and Hansen F K 2003 Phys. Rev. D 67b3001W 

Challinor A et al 2003 NewAR 47 995C 

Szapudi I et al 2001 ApJ 561 LllS 

Hivon E et al 2002 ApJ 567 211 

Balbi A et al 2002 A&A 395 417b 

Hansen F K, Gorski K M and Hivon E 2002 MNRAS 336 13041i 
Delabrouille J, Cardoso J F and Patanchon G 2003 MNRAS 346 1089 
Hinshaw G et al 2003 ApJS 148 135 
Polenta G et al 2002 ApJ 572 L27 
Komatsu E et al 2003 ApJS 148 119 
Santos M G 2003 MNRAS 341 623 

Billingsley P 1968 Convergence of Probability Measures (J.Wiley) 

Borodin A and Salminen P 1996 Handbook of Brownian Motion: Facts and Formulae 

(Birkhauser) 
Montroy T E et 2003 NewAR 47 1057 
De Gasperis G et al 2005 A&A 436 1159 



